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We study the low-frequency behavior of the ac conductivity cr(u) of a two-dimensional fermion 
gas subject to a smooth random potential (RP) or random magnetic field (RMF). We find a non- 
analytic oc \lu\ correction to Re a, which corresponds to a 1/t 2 long-time tail in the velocity correlation 
function. This contribution is induced by return processes neglected in Boltzmann transport theory. 
The prefactor of this |w|-term is positive and proportional to (d/l) 2 for RP, while it is of opposite 
sign and proportional to d/l in the weak RMF case, where / is the mean free path and d the disorder 
correlation length. This non-analytic correction also exists in the strong RMF regime, when the 
transport is of a percolating nature. The analytical results are supported and complemented by 
numerical simulations. 

PACS numbers: 73.50.Bk, 05.60.Cd 



I. INTRODUCTION 

Within the conventional approach based on the Boltz- 
mann equation, the ac conductivity of a two-dimensional 
electron gas (2DEG) is described by the Drude formula 
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where r is the transport time, v the density of states at 
the Fermi level, vp the Fermi velocity, and D the diffu- 
sion constant. Equation (|l]) corresponds to an exponen- 
tial falloff of the velocity correlation function in the time 
representation, 
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This exponential behavior of (v(t)v(O)) reflects the 
Markovian character of the Boltzmann equation descrip- 
tion and leads to analytical behavior of <7d(w) at to — > 0. 
It has been known for almost three decades, however, 
that these features result from approximations made in 
the derivation of the Boltzmann equation and are not 
generally shared by the exact solution of the problem. 
More specifically, it was shown [|) (see ^] for a review) 
that in the Lorentz gas model, where a particle is scat- 
tered by randomly located hard discs of radius a and 
density n s , there exists a "long-time tail" of the velocity 
correlation function, which has the following form in two 
dimensions (2D) in the limit n s a 2 <C 1: 
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It leads to a correction to the Drude conductivity, which 
is non-analytic at u> — > 0, 

ARecr(w) = cr o^-pM T = cr °^M T ' M T < 1 > ( 5 ) 



where I = vpT is the mean free path and we substituted 
I = 3/8n s a, the expression valid for the hard disc model. 
We will refer to this type of behavior of <j(lo) as a "clas- 
sical zero-frequency anomaly" . The long-time tail [which 
is of the form t _ ( d + 2 )/ 2 in d dimensions] can be traced 
back to processes of return of a particle to a region of 
extension ~ I around the starting point after moving dif- 
fusively during the time t ^> r || . These return processes 
give rise to non-Markovian kinetics and are neglected in 
the Boltzmann equation. 

After the discovery of weak localization the research 
interest has shifted from the above (purely classical) ef- 
fects to quantum corrections to the conductivity. For a 
non-interacting 2D Fermi gas the quantum (weak local- 
ization) correction is given by 
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where kp is the Fermi wave vector. The weak localiza- 
tion correction is of special interest, in particular, since it 
is divergent at zero frequency, indicating the crossover to 
strong localization. However, for weak disorder, kpl S> 1, 
the strong localization is of purely academic interest, 
for its observation would require exponentially small fre- 
quency and temperature and exponentially large system 
size. 

In recent years, there has been a revival of interest 
in semiclassical transport properties of 2DEG. This is 
motivated by the experimental and technological impor- 
tance of high mobility semiconductor heterostructures, 
where impurities are located in a layer separated by a 
large spacer d ~ 100 nm from the 2DEG plane. The 
(screened) random potential (RP) V(r) produced in the 
2DEG plane by the statistically distributed charged im- 
purities (density n{) is characterized by the correlation 
function HV(r — r ') = (^( r )^( r '))i which has in mo- 
mentum space the form 
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where m is the particle mass. For kpd 3> 1 (which is 
well satisfied for the high-mobility samples), the poten- 
tial varies smoothly in space and can be treated in semi- 
classical terms. Such a random potential is different from 
the Lorentz gas model in an essential way. It is weak ev- 
erywhere and shows close-to-Gaussian fluctuations (since 
at riid 2 ^> 1 potentials produced by adjacent scatterers 
strongly overlap), whereas in the Lorentz gas the poten- 
tial is zero outside and is infinite inside scatterers. There- 
fore, the Lorentz gas results cannot be directly applied 
to the 2DEG, and the problem has to be reconsidered for 
a realistic model of the random potential. 

Transport in a smoothly varying random magnetic field 
(RMF) is also of major interest. One of the main mo- 
tivations comes from the relevance of this problem to 
the composite-fermion description of a 2DEG in a strong 
magnetic field in the vicinity of half-filling of the low- 
est Landau level (y — 1/2) g. Exactly at v = 1/2 
the composite fermions move in an effective magnetic 
field B(y) with zero average and impurity-induced spa- 
tial fluctuations characterized by a correlation function 
W B {r - r') = (B(r)B(r')) of the form 



neither a clear derivation nor explicit results for the long- 
time tail in a smooth RP in the literature, we decided to 



W B (q) = (2hc/e) 2 n t e-^ d . 
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A real long-range RMF can also be realized in semi- 
conductor heterostructures by attaching superconducting 
HQ or ferromagnetic p|-pTj| overlayers or by prepattern- 
ing of the sample (randomly curving the 2DEG layer) 
PI . The strength of an RMF can be conveniently char- 
acterized fll2] , [i"3[ by a dimensionless parameter a = d/R®, 
where i?" = vpmc/ 'eBq is the cyclotron radius in a typical 
field B = (B 2 ) (magnitude of the RMF fluctuations) . 
Within the composite-fermion theory of p] this parame- 
ter is found to be equal to 1/ y/2, if the density of ionized 
impurities is assumed to be equal to the electron den- 
sity n and if correlations between the impurity positions 
are neglected. Experimental data for the magnetoresis- 
tivity around v = 1/2 are well described by the theory 
jTJJ with somewhat smaller a ~ 0.2-^0.35 (the deviation 
can be presumably attributed to the Coulomb correla- 
tions in positions of impurities and other possible features 
related to technological details of the sample preparation, 
as well as to approximations in the composite-fermion 
theory) . We will concentrate in the main part of this pa- 
per on the weak RMF case, a <C 1, when the transport 
is of conventional diffusive nature. The long-time tail in 
the case a> 1 (snake-state transport) will be discussed 
in Sec. ElIC. 



The following historical remark is in order here. After 
the initial paper jlj] by Ernst and Weyland on the Lorentz 
gas model, the i~( d + 2 )/ 2 tail in the velocity correlation 
function of a gas of particles scattered by static impuri- 
ties has been discussed in a number of publications, see 
in particular |15| , |l6| ]. However, since there appear to be 



present this material in a self-contained form (Sec. [I A 
and III A). In fact, our Eq. ( |l6| ) can be obtained from 
the mode-coupling formalism of Ref. ]l7| ]; however, the 
authors of that paper concentrate on the critical regime 
of the metal-insulator transition and do not consider the 
ac conductivity in the conducting phase explicitly. As to 
the RMF problem, which constitutes the main focus of 
the present paper, we are not aware of any treatment of 
the classical non-analytic correction to the ac conductiv- 
ity in the literature. 

For later use, we recall here the transport scattering 
rate entering the Drude formula ([!]), (||), which, in the 
case of weak long-range disorder, is found to be p8[ 



dqq 2 W v (q) (RP) , (9) 
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Let us note that the Drude result is valid in the quan- 
tum regime as well as in the classical limit. This is not, 
however, the case for corrections to the Drude result and 
we will therefore consider both the quantum theoretical 
treatment and purely classical description, for different 
parameter ranges. 

We recall that the classical description of a quantum 
particle moving in an RP or RMF characterized by a sin- 
gle spatial scale is a good approximation if two conditions 
are satisfied: (i) the quantum mechanical wavelength of 
the particle should be less than the characteristic length 
d of the disorder, i.e., kpd 3> 1; (ii) the particle should 
move incoherently, i.e., the length over which it is scat- 
tered out of its initial quantum state should be less than 
d: vft s <C d, where t s is the single-particle life time. 
The latter condition requires the random field to be suf- 
ficiently strong. We will not address the regime kpd » 1, 
vpT s 3> d in this paper. Our choice of models of disorder 
is motivated by actual physical realizations and by con- 
siderations of calculational feasibility. We will consider 
either a long-range RP or RMF governed by Gaussian 
statistics. In the quantum mechanical calculation, we will 
assume the transport scattering rate to be dominated by 
a white-noise RP. 

As discussed above in the context of the Lorentz gas, 
the fact that a particle may revisit a given region of an 
RP/RMF after large time t gives rise to a zero-frequency 
anomaly of the conductivity. Due to this effect, the ve- 
locity correlation function acquires a power-law behav- 
ior oc t~ 2 at long times, leading to Reer(tj) oc \u>\. The 
strength of this anomaly depends on the probability of 
return into the region over which the RP or RMF is cor- 
related and may thus be expected to be proportional to 
a power of d/l (the power depends on the mechanism of 
scattering) . It follows that the magnitude of the classical 
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zero-frequency anomaly for the case of smooth disorder 
considered here may be much larger than the weak local- 
ization correction (^|) in a broad frequency range. For a 
short-range potential, UfCL < 1, the role of the correlation 
length is played by the Fermi wavelength. 

We turn now to the calculation of the low-frequency 
correction to the Drude law induced by the return pro- 
cesses. 



II. RETURN PROCESSES IN THE 
QUANTUM-MECHANICAL DIAGRAM 
TECHNIQUE 




+ 




+ c.c. 



FIG. 1. Contribution to the conductivity due to return pro- 
cesses in the presence of a long-range RP. The thin dashed 
lines (forming a diffuson) correspond to a white-noise poten- 
tial, while the thick dashed line describes the long-range RP. 



We begin by considering the case in which the leading 
contribution to the transport scattering rate is given by a 
white-noise RP, while an additional weak long-range RP 
or RMF induces correlations determining the long-time 
tail. In this situation, a quantum-mechanical treatment 
of the problem is appropriate. Apart from the theoret- 
ical convenience, such a model with two types of disor- 
der is also of experimental relevance. Indeed, in essen- 
tially all realizations of real RMF (as opposed to the fic- 
titious RMF in the composite-fermion model) the trans- 
port scattering rate is dominated by a random potential 
with a relatively short correlation length (much shorter 
than that of RMF). 

A white-noise RP is characterized by the correlation 
function W(q) = (2-kut w )~ x , where 1/t w is the corre- 
sponding scattering rate (we put here H = 1). We keep 
the notation 1/r for the scattering rate related to a long- 
range RP or RMF; 1/r <C 1/t w . The diffusion process is 
described by a sum of the ladder diagrams (a diffuson) : 



r(q,w) = 
where D = v F r w /2. 



1 



1 



2niyT 2 Dq 2 — ilo 



(11) 



A. Long-range random potential 

The contribution to the conductivity induced by return 
processes is given by the sum of the diagrams shown in 
Fig- 0, yielding 

° 2 r d 2 q 



Act 
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Here S x are the vertex parts which are represented in 
Fig. @ and are given by the following expression 



S x — 
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where G^ F (p) and Gf F (p) are the retarded and advanced 
Green's functions at the Fermi energy e F . 



S x 




+ 




FIG. 2. Vertex parts of the diagrams shown in Fig. [jj The 
wavy line denotes the diffuson. 

The behavior of the correction Acr(w) at low u> is gov- 
erned by small momenta, q ~ (ui/D) 1 / 2 , in the integral 
(H). Therefore, we can make a small-g expansion of the 
vertex part (|l~3|). Expanding the integrand of ([l3|) up to 
terms linear in q pffl, we get 



S*(q) 



(14) 



Note that a naive estimate of the linear-in-q term would 
give S x (q) ~ q x ^F T w, but the two diagrams of Fig. || 
cancel each other in this order and one has to go to the 
next order in l/epT w . Substituting ([l4|), jTl|), (Q) into 
Eq. Jl^), approximating the correlation function Wv(q) 
for small q by its zero-q value, and neglecting the w- 
independent part, we find the following w-dependent con- 
tribution to the conductivity 

Aa(cj) = <7 2 - — \n(uvT w ) , \u\t w < 1 . (15) 
AePi itt 

The correction to the real part of the conductivity has 
therefore the form 



AReaH = a ^§-\u;\T w , Mr* « 1 . (16) 
oe F t w 



The condition of validity \lu\t w <C 1 given above cor- 
responds to the case d < l w . In the opposite regime, 
d^> l w , the formulas (|l5|) , (^6|) still hold, but the condi- 
tion of their validity changes to \ui\ <C D/d 2 . The same is 
valid for all the formulas for the non-analytic correction 
that are given below. 
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We first consider the situation with only the white- 
noise potential present, so that Wv{q) = WV(0). In this 
case, Eq. (rL6h yields 
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We see that the correction is small as (kplw)~ 3 , i.e., much 
smaller than the weak-localization correction (|J), and 
is therefore of minor interest. This conclusion changes, 
however, when we return to the problem with the long- 
range potential (Q) present. Equation ([lij]) then gives 
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Now the correction does not contain the quantum small 
parameter {UfIw) > which is replaced by the classical 
quantity d/l w . This prompts the expectation that the 
|cj| anomaly in cr(u>) should be essentially a classical phe- 
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nomenon. We will demonstrate it explicitly in Sec 
by calculating ARecr^) in the classical limit, where a 
long-range RP constitutes the only type of disorder in 
the system. Note that the classical limit requires two 
conditions to be met: fcpd > 1 for all relevant types 
of scatterers and also Wy(0) ^ {hvp) 2 (the latter con- 
dition means smallness of the diffraction smearing of a 
typical scattering angle; otherwise, it can be rewritten 
as vft s <C d, where t s is the single-particle life time); 
whereas Eq. (|l8| ) is obtained in the perturbative (Born) 
limit WV(0) <C (Hvf) 2 under the condition that the dif- 
fusion is due to short-range scatterers. It is also worth 
mentioning here that, in view of t > r m and vpr 3> d, 
the correction ( |l8|) is always small, ARe<7(u)/ero "C 1, in 
the range of its validity [specified below Eq. ([h])] . 



B. Long-range random magnetic field 

We consider now the same problem but with the long- 
range RP replaced by a long-range RMF. Similarly to 
( |l2| ) , we have a return- induced correction to the conduc- 
tivity 
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where 
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is the vector potential correlation function. The vertex 
part s22 + Si 2 2 is now g 
diagrams shown in Fig. ||. 



part S X a + S^xa is now given by the sum of the three 
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FIG. 3. Vertex parts of diagrams in a random magnetic 
field. 



Evaluating the vertex part at small q, we find that the 
diagrams Sx2 and Sx2 cancel each other in the order q° 
and the result is of the order of q 2 : 

(21) 



+ 42 = q 2 e F rl5 xa 
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Substituting this expression in ( |19|) and neglecting an ui- 
independent part, we find 



=- 2 MT W ) WT„«1 
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Using the explicit form of the correlation function (g), 
we get 



ARecr(cj) 7r d 
a~o ~ ~27 



\u\t w , 



(23) 



where I = vft is the mean free path characterizing the 
RMF. We see that the non-analytic conductivity correc- 
tion (and correspondingly the long-time tail of the ve- 
locity correlation function) has the opposite sign as com- 
pared to the RP case, Eq. (|lg|). This is a general feature 
of the corrections induced by a weak long-range RMF, 
as will be confirmed in Sec. IIIB by a classical calcula- 



tion for the case when such an RMF constitutes the only 
source of disorder. 



III. PURELY LONG-RANGE DISORDER: 
CLASSICAL CALCULATION OF THE 
LONG-TIME TAILS 



Having understood the nature of the \oj\ anomaly at 
the level of the quantum-mechanical diagram technique 
in the particular limit where the transport scattering rate 
is dominated by a white-noise potential, we turn to the 
case of purely long-range disorder (RP or RMF). In this 
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situation, the quantum-mechanical calculation is com- 
plicated and a classical evaluation of the non-analytic 
correction is more appropriate; it will also allow us to 
demonstrate explicitly that the correction is of classical 
origin. We will employ a formalism similar to the one 
used in for the calculation of the magnetoresistivity. 
At the quasiclassical level, the fermion gas is character- 
ized by a distribution function f(t,r,cj>), where is the 
polar angle of the velocity. The equilibrium distribution 
function is fa = 6(cf — e), where 9 is the step function. 
The deviation Sf(t, r, <p) from the equilibrium induced by 
an (infinitesimally small) external electric field Ee -4 "* 
has the form df(t,r,<j>) = eEv F ^e~ lut g{uj, r,<j>), with 
g{u>, r, (/)) obeying the Liouvillc equation 

(L + SL)g(uj, r, 4>) = cos (</> - fa) ; (24) 
Lq = —ioj + upnV . (25) 



b) 6L(r ) 6L(r') 




Here 4>e is the polar angle of the electric field and 
n = (cos <p, sin <p) the unit vector determining the ve- 
locity direction. The term Lq in the Liouville operator 
corresponds to the free motion, while SL describes the 
disorder (RP or RMF). The current density is given by 

j = — e J n^fki vSf, yielding the longitudinal conductiv- 
ity 



FIG. 4. Classical diagram technique: a) free propagator; 
b) disorder correlation function; c) diagrammatic expansion 
for the memory function M; d) first-order diagram for the 
memory function representing a return process; wavy line 
corresponds to the diffusion propagator Pd, Eq. (|3l|); e) sec- 
ond-order diagrams describing return processes, which give 
the leading contribution to the ret urn- induced correction to 
M in the random-potential case. 
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A. Long-range random potential 



Expanding ( |26| ) in SL, averaging over the RP or RMF 
(which is implicit in SL), and resumming the series, we 
get the ac conductivity in the form 



-ILJ + M 



(27) 



Here M is the self-energy (so-called memory function), 
which can be conveniently represented within a classi- 
cal diagrammatic technique (similar to the one used in 
Rcf. |16|), Fig. ||. To leading order, M is given by the 
first diagram of Fig. ^c, 



M 



-2 / — cos< 
2?r 



SL—SL > cost 
Lq 



(28) 



reproducing the results (||), ( fL0| ) for the transport scat- 
tering rate (see below) and, correspondingly, the Drude 
formula (Q). Corrections to the memory function M(oS), 
which correspond to the return processes, are evaluated 
in Sec. |IIIA| and |IIIB| for the cases of RP and RMF, 
respectively. 



The fluctuating contribution to the Liouville operator 
due to the RP is found to be 



SL V = <5v(r)nV + (V£u(r))nj 







(29) 



where nj^ = z x n = (— sin tfi, cos <f>) and Sv (r) = v (r) — vf 
is the deviation of the local velocity v(r) = {(2/m)[eF — 
y(r)]} 1 / 2 from its average value vf- The leading-order 
contribution ( p8|) to the memory function reads 
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reproducing the transport scattering rate defined by 
Eq. (||), Mo = 1/r. The first-order diagram describing 
the return process is represented in Fig. |]d. The cor- 
responding expression is obtained by replacing the free 
propagator 1/Lq in (|30|) by the diffusion propagator 
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7(q, <^)7(q 5 <^') 



Dq 2 



IOJ 



(31) 



1 — iql cos(0 — (f) q ) , ql -C 1 



5 



The replacement yields the return-induced first-order 
correction to the memory function 



d 2 q 



sin m sin ( 
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x q 2 W v (q)P D (( l , <j>, <t>') smcf) 1 sin {dp 1 - fa) . (32) 

Evaluating the u> dependent part of ( |32| ) at ojt <C 1 and 
approximating (as in the quantum-mechanical calcula- 
tion) Wv(q) by its value at q = 0, we find 



AMi(w) _ HV(0)ut 



Ml 



which gives 
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(34) 



for the specific form (]?]) of the correlator W*V- 

Now let us show that, in actual fact, the leading con- 
tribution to the non-Markovian correction to ReM(w) 
comes from second-order processes described by the two 
diagrams in Fig. ||e, whereas that given by Eq. ( ^3|) can 
be neglected in the first approximation. Specifically, the 
second-order term AReM2/Mn ~ (d/l) 2 \uj\T scales with 
a smaller, as compared to AMi, power of the param- 
eter d/l -C 1, despite having one more impurity line. 
This, at first glance, counterintuitive feature is related to 
the anomalous smallness of AMi in the otherwise "regu- 
lar" expansion in powers of d/l (third- and higher-order 
terms in AM can be shown to be negligible compared 
to AM2). We first explain this feature by using the fol- 
lowing power-counting argument. The |w|-anomaly in 
AM comes from the integration over small q of the form 
J d 2 q q 2 J (Dq 2 —iuj), where the numerator of the integrand 
tends to zero as q 2 at q — ► 0. In Eq. d32|), the factor of 
q 2 is related to the vanishing of the correlator 



d r exp(— iqr) (VSv(0)VSv(r)) oc q 2 W v (q) 



in the limit q — ► 0, since the correlator carries the small 
momentum q, the same as the diffuson, according to 
Fig. |]d. On the other hand, in second order in Wy, the 
large momenta flowing through impurity lines are "dis- 
entangled" from the small momentum q carried by the 
diffuson. The leading q 2 -term comes now from the 0(ql) 
corrections to the diffusion propagator given by the fac- 
tors 7(q, 4>) in Eq. (|3l]). Let us count powers of I: two 
factors 7(q, cf>) yield q 2 l 2 , whereas one loses only I -1 when 
going to the second order, which explains the total gain 
of one power of l/d as compared to Eq. (34). 

The expression for AM2 at u> — > obtained from the 
sum of the two diagrams in Fig. [|e reads [we neglect the 
dependence on q everywhere but in Po(q, <t>, </>')] 



4i fdd>dd>' d 2 q d 2 k , ~ „ , s 
x cos <j)A(k, 4>)P D (q, <f), <p")Im A(k, 0') cos 0' , (35) 



where 
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We thus obtain 
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Since 



ARea(w) AReM(w) _ ARep(w) 
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where p{oS) — a 1 (w) is the ac resistivity, we get finally, 
using Eq. @) for W v , 

ARea(w) 3tt f d\ 2 , , 

_U = T ^ 7 j Mr, Hr«!. (39) 

The prefactor of the |a; [-correction to Recr(a;) is posi- 
tive, as in the quantum-mechanical result (|l^) and in 
the Lorentz gas formula (||). Note that the correction 
(|39l ) matches that for the Lorentz gas (||) at n s d 2 ~ 1, 
as expected — since this condition separates two extremes 
of strongly non-Gaussian (Lorentz gas) and Gaussian 
[Eq. (Q)] disorder. On the other hand, the crossover be- 
tween Eqs. (H) and © occurs when the two following 
conditions are fulfilled: kpd ~ 1 and W-V(O) ~ (frvp) 2 [cf. 
the definition of the classical limit for Gaussian disorder 
after Eq. ©]. 



B. Long-range random magnetic field 

The fluctuating contribution to the Liouville operator 
induced by the RMF has the form 



5L B = —B{v)^- 

mc 00 



(40) 



The lowest-order contribution (|28|) to the memory func- 
tion 
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reproduces again the corresponding transport scattering 
rate (]l0|), Mo = 1/t. The first-order correction due to 
return processes, Fig. |]d, reads 
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2 (-) 

Vmc/ 



d(f>d(p' d 2 q 



sin < 



2tt 2tt (2tt) 2 
xW B {q)P D {q,q>,q>')sm<f>' , 



(42) 



where the factors 7(q, (j>) should be included in 
Po(q, cf), (/)'), which gives 



ARecr(cj) AReM(w) 



00 



Mo 



-) 



Wb(0) , 
8t$ 1 



(43) 



Note that, in contrast to the case of RP, the leading con- 
tribution to the return-induced correction AM in RMF 
comes from the first-order processes. This is because 
the RMF scattering operator SLb (in contrast to its RP 
counterpart SLy) does not involve spatial gradients. Us- 
ing the RMF correlation function (pi), we finally get 



ARecr(w) 



(TO 



7T d 

27 



lo\t 



2 \lo\t , |w|t«1. (44) 



We have found, therefore, in agreement with the 
quantum- mechanical result (p3|), a negative sign of the 
|w|-contribution to the conductivity. Analyzing the cal- 
culation, one can trace the difference in sign [as compared 
to the Lorentz gas result (||) and the RP results (p"8|), ( |39|) ] 
back to the fact that the RMF scattering operator ( |40| ) 
is odd with respect to time reversal. 



1.03 




FIG. 5. Real part of the ac resistivity in the RMF with 
a — 0.5 normalized to its lu = value. The dashed line is 
a guide for the eye, ARep(w)/p(0) cx |cj|. The positive pref- 
actor of the |tj [-correction to Rep corresponds to a negative 
prefactor for Reo\ 

To check the above analytic findings, we have per- 
formed numerical simulations of the classical motion of a 
particle in an RMF. The results obtained for the mem- 
ory function [BU at a = 0.5 are shown in Fig. [| We find 



a positive |w [-correction to the real part of the memory 
function (or, equivalently, resistivity), which corresponds 
to a negative correction to Reer, in agreement with the 
theoretical result (|44|). The magnitude of the correction 
is, however, considerably smaller than Eq. ( p4| ) would 
predict. We attribute this discrepancy to the fact that 
Eq. ([l4|) was derived for a <C 1 and, apparently, the 
numerical value of the coefficient in this formula cannot 
be trusted for a as large as 0.5 Unfortunately, at 

smaller values of a < 0.2, the effect becomes so weak that 
it is swamped by the statistical noise. A smaller value 
of the coefficient at a = 0.5 [as compared to the a <C 1 
formula (|^|)] is further consistent with the fact that at 
a > 1 the coefficient changes sign and the correction to 
the conductivity becomes positive (see below). Let us 
also note that the range of validity of the \lo [-correction 
found numerically is in full agreement with the theoreti- 
cal expectation (lot < 1). Indeed, as is seen in Fig. |[ the 
linear increase of Re p(to) holds up to lot ~ 0.65 where it 
transforms (rather abruptly) into a falloff (related to the 
ballistic motion on time scales t <t). 



C. Strong random magnetic Held: long-time tail in 
transport on a percolating network 

In a strong RMF (a » 1) the character of the trans- 
port changes drastically. In this regime, the diffusion 
takes place in a restricted space and is determined by 
a small fraction of trajectories — so-called "snake states" 
p3| , p4| — which wind around the B(r) = contours. Since 
the snake states can go over from one B(r) — line to 
another at saddle points of the RMF (where the two con- 
tours come sufficiently close to each other), they propa- 
gate effectively on a percolating network |2|,|l^,|l^] for 
which such saddle points serve as nodes. This network 
is characterized [[f4j by a typical length of a link, L s ~ 
da 14 / 9 , and a typical distance between two neighboring 
saddle-points (size of an elementary cell), (, s ~ da 8 / 9 . 
The different scaling of L s and £ s with a is due to the 
fact that the structure of the links of the network is frac- 
tal. The network is chiral, i.e., the links are directed; each 
node has two incoming and two outgoing links. Since the 
snake-state velocity is of the order of the Fermi velocity, 
a characteristic time of traversal of a link is t s ~ L s /vf- 
The quasiclassical dc conductivity in this regime was cal- 
culated in |p^| , |l4| , the result being a ~ kpd/a 1 ^ 2 In 1 / 4 a. 

In Ref. |l4| we argued, on phenomenological grounds, 
that for such a percolation-type transport problem there 
should be a non-analytic contribution to the ac conduc- 
tivity of the form 



ARe a(uo) 
"(0) 



< 1 



(45) 



Below we demonstrate how this result comes about in a 
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network model due to fluctuations in the geometry of the 
network. 

Let us start with a regular square network (Fig. ^a) 
with all links characterized by the same distance £ s be- 
tween the end points and by the same "flight time" r s . 
We also assume that the probability of turning in either 
of two allowed directions at each node is 1/2. The classi- 
cal diffusion constant is then D = £^/At s . It is straight- 
forward to see that there is no memory effect in classical 
transport on the regular network: the velocity correlation 
function is exactly zero for t > t s . Let us now study the 
effect of fluctuations in the network geometry, i.e., in vec- 
tors connecting the beginning and the end of individual 
links. Since to describe such fluctuations quantitatively 
in a real percolating network is hardly possible, we con- 
sider the following model. We imagine the regular lattice 
considered above perturbed by a small fraction <C 1 of 
"defects" of the type shown in Fig. |6|b (a defect can have 
any of four possible directions). We assume the flight 
times of all links to be equal (we will discuss the effect of 
fluctuations in the flight times later). 



a) 



(v(nr.)v(O)) = —J2 E E v^v^an-ljr.) 



AN ^ ^ ^ 

ij ji=l,2 w=3,4 



(46) 



where N is the normalization factor (total number of 
sites), Pi jit) is the probability of moving from a site i to 
a site j in a time t, and = ^^/t s is the velocity at 
the link (ifi). The majority of the sites i,j will give zero 
contribution to ( |46"| ) after the summation over fi and v, 
since the velocities of the two outgoing (or two incoming) 
links are exactly opposite to each other for the regular 
lattice. A non-trivial contribution will come from terms 
with both i and j lying at a defect. Indeed, consider the 
term with i = 1, j = 2 (Fig. ^]b). The corresponding 
contribution to (Ehf) is 



n d Pi2{t - T s ) 



(47) 



The probability density in a continuum model for a dif- 
fusing particle to move a distance r in a time t is 



P(t,r) = 



1 



2 /4Dt 



AirDt 

Therefore, the probability Pv2,[t) for t ^> t s is 

"T, 



Pl2(t) 



irt 



1 + 



(48) 



(49) 



This return process yields a contribution to the velocity 
correlation function of the form 



tsY Ts 1 



b) 




FIG. 6. Chiral network model: a) regular network; b) de- 
fect on the lattice. 

For each lattice site j, we label adjacent links as (j/J,), 
with fj, = 1,2 for incoming and fj, = 3, 4 for outgoing 
links. The velocity-velocity correlation function for a 
time t — nr s (with an integer n) can be written as 



However, this 1/i-contribution is canceled if we take into 
account the terms with i = 3, 4 and j = 5, 6 as well. The 
total contribution reads 

(v(i)v(O)) = n d [P 12 {t - t s ) - P 15 (t - t s ) 

-P 32 (t - t s ) + ip 35 (< - t s ) + ^P 36 (t ~ t s )} . (50) 

Since all of the relevant return probabilities have the 
form jigj), the l/t terms cancel. It is easy to see that 
this cancellation has a general character, i.e., is indepen- 
dent of the particular structure of the defect. We thus 
conclude that the result is of the next order in r s /i, 

(v(t)v(0))~-n d (|) 2 (^) 2 . (51) 

While we do not calculate the numerical coefficient in 
( pl| ) pq |, we see no reason which would require it to be 
zero, so that we believe that it is generically non-zero. 
Setting now rid ~ 1 for a realistic (strongly fluctuating) 
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network results in a non-analytic correction to the con- 
ductivity of the form (^). As to the sign of the effect, 
we have to resort to numerical simulations (see below). 

Let us now consider the effect of fluctuations in flight 
time. We return to the regular square lattice (with the 
lattice constant £ s ), but now allow for variation of the 
flight times t m from one link to another. We will show 
that in this model the 1/t 2 tail does not exist. Equation 
( [4(f ) for the velocity-velocity correlation function is now 
modified as follows 



v(*Mo)) = ^EE E 

s ij /i=l,2 v=3,4 
dT , Wif l Wj v Pij{t -t-t'] 



{ I dr 

>0 -h\ 



(52) 



where = £^/t^ and r s = (r M ). Since the fluctuations 
of the flight times of different links are uncorrelated, a 
non-zero contribution to ( |52"| ) comes only from neighbor- 
ing sites i,j connected by a link going from j to i [in 
other words, one of the links (z/u) should be identical to 
one of the links (jv)\- We find thus 

(v(i)v(O)) - 

dr / dr' P lj (t-T-T r )) 




(i 



tit 2 



dr [ dr 1 '\Pij{t - t - t') 
Jo T i / Tl 



(53) 



Expanding Pij(t — r — r') ~ 1 /7r(t — r — r') in r and r', 
we see that the terms of the 1/t and 1/t 2 orders cancel, 
and the leading non- vanishing contribution is of the 1 /i 3 
order, so that the corresponding contribution to Re cr(w) 
shows a weak non-analyticity oc to 2 In \u>\ only. Note that 
for a non-directed network, fluctuations of the flight time 
yield a still weaker non-analyticity ARe <j{uj) oc |w| 3 , or, 
equivalently, a 1/t 4 long-time tail |37fl . 

Since a real percolating network exhibits all possible 
sorts of fluctuations, the fact that we find the 1/t 2 tail in 
the model with fluctuating £ M 's is sufficient to conclude 
that such a tail should be present in the problem of the 
transport in strong RMF. In Fig. [?] we show the results 
of the numerical simulations of the problem for a ~ 4. A 
pronounced dip in the ac conductivity around uj = in 
the expected range of frequencies \ui\ < l/r s ~ vf / 'da 14 / 9 
nicely confirms our analytical conclusions. The sign of 
the non-analytic correction corresponds to a decrease of 
Recr as |cj| — * 0. 

It is worth mentioning that the problem of a random 
walk on such a percolating network is a close relative 
of the advection-diffusion problem in a spatially random 
velocity field v(r) ("steady flow") with V-v = ("incom- 
pressible liquid") characterized by the correlation func- 
tion 



d 2 r exp(-z'qr) (v a {0)v/3(r)} = W v (q)(S al3 q 2 - q a q/3) 



(54) 



This model was studied in a scries of papers |28| |3C|] 
with emphasis on the case of long-range correlations, 
namely W v (q) oc q~ 2 for q — > [which corresponds 
to (v a (0)vp(r)) oc r~ 2 ]. In contrast, we have consid- 
ered a percolation lattice with short-scale distortions 
[W v (q) — > const at q — > 0] . One can check (see Appendix) 
that the advection-diffusion problem yields a t~ 2 tail in 
this case, in agreement with our consideration above. 
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FIG. 7. Real part of the ac conductivity in the RMF at 
a = 4.04. A non-analytic dip around u> = is clearly seen. 
The low-frequency increase of conductivity is restricted to the 
region \lo\L s /vf < 1 (where L 8 ~ da 1A ^ is the length of a link 
of the percolating network), in agreement with the theory. 



IV. EFFECT OF INELASTIC SCATTERING 

So far our considerations have not included inelastic 
scattering processes which change the energy of a par- 
ticle. The question arises whether the zero-frequency 
anomaly ARe a(ui) is cut off at low frequencies w ~ 1 /T- m , 
where r; n is a relaxation time for the inelastic processes. 
This question is studied most conveniently within the 
Liouvillc-Boltzmann approach of Sec. [II. To this end, we 
consider the linearized distribution function 5f(w, r, e, <fi) 
of particles with energy e and velocity direction specified 
by the polar angle <j>, subject to a smooth RP or RMF 
and inelastic collision processes, obeying the Liouvillc- 
Boltzmann equation 



{-iuj + vV + 6L)6f - I in {Sf) = S 



(55) 



with the source term S = evE(<9 /q / de) . Here /o is the 
Fermi distribution function corresponding to a tempera- 
ture T (which we will assume to be low, T -C Ep). A 
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simple model form of the collision integral I- m , which re- 
spects particle number conservation, is sufficient for our 
purposes: 



i 



Sf(e,4>) 



dfo 
de 



dc' 



d(j)' 
2^ 



5f(e',<f>') 



(56) 



(For simplicity we adopt the model of isotropic, energy 
independent inelastic scattering.) 
The conductivity is obtained as 



2 2 

e vv F 



dc 



dfo 
' de 



, COS d>^- COS ( 

2tt \ Y Lo 



(57) 



where Lo = — iuj + Tj~ + vpiiS7 . Expanding in SL, av- 
eraging over the long-range disorder and resumming the 
series, one finds u{oj) in the form J27j). The memory func- 
tion M is now given in lowest order for the cases of RP 
and RMF by ( |30| ) and ( fll| ) , respectively, with ui replaced 

In order to calculate the effect of return processes, one 
needs to know the diffusion propagator Po(q; e, cb; e' , cb 1 ) 
for particles starting with energy e and velocity an- 
gle cb and returning with energy e' and angle cb' . It 
obeys Eq. (|55| ) with the source term replaced by So = 
5(cb—ct>')5(e — e')(—dfo/de). After averaging over 5L, one 
finds 



P D (q;e,cb;e',cb') = 



dfo\ 7 e (q,</>h £ '(q, <j>' 



de J D t q 2 
regular terms , 



dfo 
de' 



(58) 



where D t = v F r t /2 with r t = r _1 is the total dif- 

fusion constant including elastic and inelastic scattering 
processes, and 7 e (q, <p) — [1 — iqv(e)rt cos(cb — cb q )} with 
v(e) = {2e/ra) 1 / 2 . As expected, the diffusion propagator 
shows a diffusion pole even in the presence of inelastic 
processes, due to particle number conservation. The scat- 
tering "out" of particles with given energy e into other 
energy states is exactly compensated by a corresponding 
scattering-in contribution. 

Let us define the function AM(e, e') in the same way 
as in Eqs. (f35|) , (fl2]) with the only change Po(q, <j>, cb') — > 
Pc(q; e, 4>\ e ' ) 4>')- The correction to the memory function 
due to return processes, AM, is then given by 



AM = J dede'AM(e, e') . 



(59) 



As a result, the expressions (37), ( |3S| ) and (|43|), ( f44| ) 
remain valid, provided (i) r and / are replaced by the 
full momentum relaxation time Tt and mean free path 
It = VFTt, respectively, and (ii) additional factors of 
(^t/r) 2 and r t /r are included in (^), (f44|), respectively, 



which stem from the explicit factors of W in the defini- 
tion of AM(e,e'). 

Thus, the classical zero- frequency anomaly is not cut 
off at finite temperature. This should be contrasted 
with the quantum zero-frequency anomaly induced by 
the weak-localization and Altshuler-Aronov (interplay of 
interaction and disorder) effects. It follows that increas- 
ing temperature favors the experimental observation of 
the classical anomaly. 



V. CONCLUSIONS 



In this paper, we have studied memory effects in the 
low-frequency ac conductivity of a 2D fermion gas in a 
long-range random potential or random magnetic field. 
We have calculated the long-time tail in the velocity cor- 
relation function induced by diffusive returns of a particle 
and leading to a non-analytic |cj|-beIiavior of the real part 
of the conductivity (zero- frequency anomaly). While in 
a random potential the |w [-contribution is positive (as in 
the Lorentz gas, Ref. |ij) and is proportional to (d/l) 2 , 
where <g; 1 is the ratio of the correlation length to 
the mean free path, a smooth weak RMF induces a much 
larger (oc d/l) correction of opposite sign. The sign dif- 
ference can be traced back to the RMF scattering being 
odd with respect to time reversal. 

We have also demonstrated how an |w|-contribution to 
Recr(w) arises in the regime of strong random magnetic 
field, where the transport is determined by percolation 
of the snake states. In this case, spatial fluctuations in 
the geometry of the percolating network are responsible 
for the memory effects. 

Our numerical simulations confirm the existence of 
these non-analytic contributions at low frequency, as well 
as the unconventional sign of the correction in the weak 
random magnetic field. With increasing strength of the 
RMF, when the system crosses over into the regime of 
the percolating transport, the sign of the effect changes. 

The experimental observation of the non-analytic low- 
frequency behavior of the ac conductivity of composite 
fermions would be of considerable interest. In particu- 
lar, we predict that the |w|-term in Recr is negative at 
v = 1/2 in the high- mobility samples (where the strength 
of the effective RMF is a - 0.3 @), but should change 
sign if the system is driven toward the percolation regime 
by adding more long-range scatterers (e.g., antidots pl|). 
A sign change is also expected with increasing effective 
magnetic field (moving away from half-filling) , when the 
system crosses over to a percolation regime, as discussed 
in Ref. |TI. 
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APPENDIX 

Consider a particle moving in a diffusive medium with 
a diffusion coefficient D, subject to a spatially random 
velocity field v(r). Let the velocity field be incompress- 
ible (V • v = 0) and determined by the correlator ( |5^ ) 
with a finite W v (0). Assuming the random field to be 
weak, we can expand the Green's function 

G(r,r') = (r|(-wj- £>V 2 - V • v) _1 |r') 
in v(r), which yields for the Fourier transform 
1 



G(q) 



-iu> + Dq 2 



1 



(-iLU + Dq 2 ) 2 



(vV)- 



1 



-iw - DV 2 



(v-V) 



Using the correlation function (|54|) and resumming the 
series, we thus find 



G{q) = 



1 



-iw + (D + SD)q 2 ' 
with the following correction to the diffusion coefficient 



6D 



d 2 q q 2 W v {q) 
(2tt) 2 -iuj + Dq 2 



3 W v (0) 
16 D 2 1 



Hence, this continuous model predicts a t~ 2 tail in the 
velocity-velocity correlation function, with a positive co- 
efficient. Note, however, that in a lattice model the sign 
depends on the microscopic structure of disorder. 
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